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We investigate the quasiparticle density of states in disordered d-wave superconductors. By 
constructing a quantum map describing the quasiparticle dynamics in such a medium, we explore 
deviations of the density of states from its universal form (oc E), and show that additional low-energy 
quasiparticle states exist provided (i) the range of the impurity potential is much larger than the 
Fermi wavelength [allowing to use recently developed semiclassical methods]; (ii) classical trajectories 
exist along which the pair-potential changes sign; and (iii) the diffractive scattering length is longer 
than the superconducting coherence length. In the classically chaotic regime, universal random 
matrix theory behavior is restored by quantum dynamical diffraction which shifts the low energy 
states away from zero energy, and the quasiparticle density of states exhibits a linear pseudogap 
below an energy threshold E* < Ao. 



PACS numbers: 71.23.-k, 74.72.-h, 05.45.Mt 

In recent years, considerable attention has been fo- 
cused on the low-energy properties of the quasiparticle 
spectrum of disordered cuprate superconductors 0, 0- 
Because many of the cuprate superconductors are ran- 
domly chemically doped insulators and disorder is a pair- 
breaker for d-wave superconductors, the role of nonmag- 
netic impurities is particularly important for an under- 
standing of the d-wave superconducting state, its quasi- 
particle spectral and transport properties. Of special in- 
terest is the low-energy behavior of the single-particle 
Density of States (DoS) p(E). 

In early work jj] , the self consistent T-matrix approx- 
imation was shown to break down for 2-dimensional d- 
wave superconductors. This led to a series of papers us- 
ing nonperturbative methods, which predicted (at first 
sight) contradictory results: vanishing [j| 0, El , con- 
stant HQ, and diverging DoS as E -> 0. On the 
numerical side, several investigations also predicted both 
vanishing 0, 0] , and diverging 0, Q] DoS at zero en- 
ergy. It was soon argued, based on numerical analysis, 
that the reason behind these contradicting predictions 
is the fact that the microscopic details of disorder (i.e. 
details beyond the transport mean free path I, such as 
the density of scatterers or the correlation length ( of 
the impurity potential) as well as the symmetries of the 
clean Hamiltonian matter both qualitatively and quan- 
titatively 0,^1 (see also [lif). An important feature 
shared by the numerics of Refs. El El El El El is that 
the disorder is introduced via isolated point-like scatter- 
ers. Long-wavelength disorder, which may arise due to 
chemical doping away from Cu02 planes, or can be in- 
duced via ion radiation techniques 17] or via a STM 
tip E3 > is thus ignored. Effects of long- wavelength dis- 
order are expected to become dominant when the Cu02 
planes have almost no atomic disorder 0]; they are the 
main focus of this paper. 

It has recently been realized, in the context of meso- 



scopic physics and weak localization, that £ and £, to- 
gether with the Fermi wavelength Af, define two classes 
of complex quantum systems: quantum disordered sys- 
tems where Xf£/C 2 > 1 an d quantum chaotic systems 
for which Xpl/C, 2 < 1 j^. The latter class is character- 
ized by the emergence of a new diffractive scattering time 
scale, te = v~ l lnfC/Ai?], defined as the time it takes for 
the classically chaotic dynamics (with Lyapunov expo- 
nent v) to stretch a wavepacket of minimal initial exten- 
sion Xp to a length £. In contrast to quantum disordered 
systems, quantum chaotic systems exhibit nonuniversal 
properties due to their short-time classical (i.e. deter- 
ministic) dynamics. In particular, significant deviations 
from Random Matrix Theory (RMT) emerge, as was re- 
cently found by Adagideli et al. [lfj in the context of 
impurities in d-wave superconductors. These authors 
used a semiclassical approach to calculate the low en- 
ergy DoS for a collection of extended scatterers (a quan- 
tum chaotic system) and found an asymptotic behavior 
p{E) ~ 1/E\ lni? 2 | 3 as E — ► 0. They nevertheless argued 
that the RMT predictions of a linear pseudogap would be 
restored at lower energy E < E* , i.e. that the singularity 
in the DoS would be cut off at an energy E* related to te, 
by diffractive (nonclassical) scattering occurring at larger 
times r > Te- The purpose of the present paper is to in- 
vestigate the modifications that the DoS undergoes as the 
correlation length of the impurity potential increases and 
te becomes relevant. We will focus our attention on (i) 
providing for numerical checks of the theory of Adag ideli 
et al. in the case of long- wavelength disorder [lfj; (ii) 
finding out whether for some E* the DoS is suppressed 
for E < E* , in agreement with RMT predictions 0,0,0], 
(this could also reconcile the above mentioned contra- 
dictory predictions); and (iii) investigate the transition 
region between extended and pointlike disorder. 

We start by introducing a quantum map model for 
quasiparticle states in disordered d-wave superconduc- 
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tors. The main motivation behind this model is to inves- 
tigate discrepancies (in low-energy DoS) between point- 
like vs. extended disorder as well as the transition region, 
i.e. the regime in which the impurity size is intermedi- 
ate. To the best of our knowledge no disorder model 
which includes extended impurities as well as pointlike 
in d-wave superconductors has been studied numerically 
so far. This map has two additional advantages: First, 
as both the density and correlation length of impurities 
can be tuned independently, it is possible to interpolate 
between the two extreme regimes of strong disorder: uni- 
tary disorder (i.e. disorder due to dilute, pointlike scat- 
terers I El EL quantum disorder) and quasi- 
classical disorder [lfj (i.e. disorder due to extended scat- 
terers, viz. quantum chaotic). Second, from a numeri- 
cal point of view, it allows for the investigation of very 
large system sizes, i.e. lattice sizes of up to 256 x 256, 
which are necessary for both variations of the disorder 
correlation length and the numerical extraction of the 
parametric behavior of the DoS. Our reasons why a dy- 
namical model is relevant are: (i) In absence of supercon- 
ductivity, many properties of quasiparticles in disordered 
media (such as Anderson localization) are correctly de- 
scribed by 1-D maps j^J. I n fact it has been shown by 
Altland and Zirnbauer that one of those maps, the 1-D 
kicked rotator, and quasi-one-dimensional metallic wires 
are described by the same effective field theory j^l • Re- 
cent numerical investigations in D = 2 suggest that this 
is also true in higher dimensions |23| . (ii) In presence 
of superconductivity, Andreev maps based on the kicked 
rotator have recently been shown to adequately describe 
quantum dots in contact with a superconductor |24j |. 

We first briefly discuss generic properties of quantum 
maps for uncoupled quasiparticles. The dynamics corre- 
sponds to a succession of free propagations, interrupted 
by sudden kicks of period To, i.e. instantaneous perturba- 
tions. Quantum maps are conveniently represented by a 
unitary, Floquet operator F, giving the time-evolution af- 
ter p kicks as u(p) — F p u(0), for an initial wavefunction 
u(0). The matrix F has eigenvalues exp(— i£ m ), which 
define quasi-energies e m € (— 7r,7r) (energies and quasi- 
energies are expressed in units of H/to). While the en- 
ergy is not conserved, the periodicity of the kick still 
preserves quasi-energies, much in the same way as a pe- 
riodic potential breaks translational symmetry, but still 
preserves quasi-momentum. Time evolution of hole ex- 
citations (being the time-reversed of electronic excita- 
tions) is given by v(p) = (F*) p v(0). Specializing to the 
D-dimensional kicked rotator, we write the Floquet op- 
erator as |25j 
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It describes the free motion of a particle with dimension- 
less coordinates {rj} (e.g. expressed in units of a lattice 
constant), which is interrupted at periodic time intervals 



To by a kick of strength K-H?^ cos rj . I is the moment of 
inertia of the particle, and K is the kicking strength. For 
D = 1 and 2, increasing K makes the classical dynamics 
evolve from integrable (K = 0) to fully chaotic [K > 7, 
with Lyapunov exponent A w ln(K/2)]. For < K < 7 
stable and unstable motion coexist (a so-called mixed 
phase space) [25| . Increasing K is thus tantamount to in- 
creasing the amount of disorder, the fully chaotic regime 
corresponding to a finite density of impurities. 

Electron and hole excitations inside a superconductor 
are however coupled by a nonvanishing pair-potential. 
Accordingly we extend the kicked rotator of Eq. (JTJ to 
a Bogoliubov-de Gennes form. We discuss this construc- 
tion for the case D = 2. First, we replace the free quasi- 
particlc motion by a coupled electron and hole dynamics, 



To = exp(— iTi.ro/ H), 
Ti = H a z + A a x . 



(2a) 
(2b) 



Here, H = — (h 2 V 2 /2I) — E F , with E F the Fermi energy, 
<Tx,z are Pauli matrices acting in particle-hole space, and 
A is the superconducting pair potential. Second, the 
coupled quasiparticle motion is followed by a kick 



•F k = exp(-iTi K /h), 



Hk = 
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■ cos x cos y <j z . 



(3a) 
(3b) 



Exponentiating the Pauli matrices, we end up with the 
Bogoliubov-de Gennes-Floquet (BdGF) operator 



T — Tk Fo, 



To = cos v/(i/ 2 + A2)(r A)2 J 
isin v / (H 2 + A 2 ){r /h) 2 
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cos x cos y 



(4a) 

Hcr z + A(j x }, (4b) 
(4c) 



with X, the identity matrix in particle-hole space. For 
A = 0, Eq. (0} describes uncoupled electron and hole 
excitations in a disordered 2D metal. Once this metal 
becomes superconducting, A couples these excitations 
during their free propagation, while it is neglected dur- 
ing the instantaneous kick. As in the case of a BdG 
eigenproblem, the 2M quasienergies [with average spac- 
ing S = (e m +i — £m) — ir/M] of the BdGF equation 
T <f> m = exp(—ie m )(pm, come in pairs with opposite sign 
e m = — £2M— m+ii similarly to the spectral properties of 
a BdG Hamiltonian. These considerations establish the 
correspondence between the map of Eq. and quasi- 
particles in a dirty superconductor. 

We next quantize the phase space on a 4-torus 
{x, y,Px,Py}, with dimensionless momentum p xy = 
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FIG. 1: Density of states for the d-wave (right) and extended 
s-wave (left) Bogoliubov-de Gennes kicked rotator defined by 
Eqs. 12IH . and parameters L x x L y = 128 x 128, Ao = 0.4, 
Ef = 27r 2 /5, and K = 0. (solid lines), 2. (dashed lines) and 
8. (dotted-dashed lines). 

-ih e $d/d(x,y) £ (0,2tt) The effective Planck con- 
stant h e g = Htq/Iq takes on values h e g = 2ir/M, with 
integer M — L x ■ L yi in term of the real-space linear 
system sizes L x ,y (also expressed in units of a lattice 
spacing), and the impurities have a spatial extension 
( = 0(L x ,L y ). The BdGF operator is then a 2M x 2M 
unitary matrix, and we consider the two cases of d-wave 
[A(p x ,p y ) = A (pI-pI)/(pI+pI)} and extended s-wave 
[A(p x ,p y ) = A \p 2 x ~p 2 y \/(pl + p 2 )] pair potentials, for 
which J-q is diagonal in momentum representation. Not- 
ing that Tk is diagonal in real space representation, we 
rewrite T as 



j ppf 



(5) 



where U = UT, and U is the unitary matrix of the 2D 
Fourier transform between real space and momentum co- 
ordinates, Upffi = M -1 ' 2 exp[(27ri/M) p-p'}. We numer- 
ically extract the quasienergy DoS from the eigenvalues 
shi£ m of the hermitean matrix -^{T — which we 
diagonalize using the Lanczos algorithm 26] . 

In Fig. 1 we show the quasienergy DoS for d-wave and 
extended s-wave pair potentials away from half filling 
(Ef = 2ir 2 /5 < 7r 2 /2), as the kicking strength increases. 
In the clean case (K — 0) the two DoS are the same. 
The gap singularity at E/Aq — 1 gets washed out as K 
increases in both cases, however, a peak emerges in the d- 
wave DoS around E = 0, while p(E) = in the extended 
s-wave case. This is in agreement with Ref. 0,H3, i.e. 
the existence of low-energy states requires a change in 
the sign of the pair potential. In the extended s-wave 
case, the low energy peak is shifted by an energy corre- 
sponding to the gap averaged over all momenta mixed by 
the impurity potential. 

We focus on the d-wave symmetry from now on. A 




FIG. 2: Main plot: Low energy DoS for the d-wave 
Bogoliubov-de Gennes kicked rotator of Eqs. J21 1 1 . 
L x x L y = 256 x 256, A = 0.4, E F = 2tt 2 /5, and 
K = 8. Inset: Asymptotic of the density of states at 
low excitation energy for the same set of parameters, for 
long-wavelength disorder (quantum chaotic; black circles) 
and diffractive disorder as defined in Eq. @ with Nh = 27 
(quantum disorder; empty diamonds). The solid and dashed 
lines indicate a oc E and oc _E _1 behavior respectively. 

closer look at the DoS in the fully chaotic regime with 
K = 8. is provided in Fig. 2. It indicates that the char- 
acteristic semiclassical singularity exhibited by the DoS 
as E — * is cut off at an energy £* < Ao, where a sharp 
drop occurs and p(E) — > 0. RMT predicts such a drop 
to occur over an energy scale given by the Thouless en- 
ergy 0], which in our case is however significantly larger 
than Ao jl^]. We thus attribute this drop to the emer- 
gence of diffractive scattering at times larger than te 
as follows. According to Ref. 0, the DoS correspond- 
ing to low energy semiclassical states can be estimated 
from a mapping onto a tight-binding chain with random 
hoppings, for which the eigenfunctions are localized with 
an energy-dependent localization length (;(E) [2^. At 
low energies, £ exceeds the diffractive scattering length 
vfte, (vf is the Fermi velocity) in which case hoppings 
between otherwise uncoupled tight-binding chains (cor- 
responding to different classical trajectories) have to be 
taken into account. The emergence of these processes sig- 
nals the breakdown of semiclassics and the restoration of 
RMT. One thus expects the vanishing of the DoS below 
a threshold energy given by the condition £(E*) ~ vfte- 
Since £ (E) is bounded b y th e superconducting coherence 
length, £(E) > Hvf/A 10], the observation of the semi- 
classical peak in the DoS requires a long enough diffrac- 
tive scattering length vfte > ^w_f/A . While prelimi- 
nary results corroborate this argument, a detailed inves- 
tigation of E* will be presented elsewhere |3fT |. 

In the inset to Fig. 2 we show the asymptotic behav- 
ior of the DoS on a log-log scale. Once abstraction is 
made of the drop in the DoS below E* , the semiclassical 
data exhibit a singular behavior slightly below oc E^ 1 
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FIG. 3: Low energy DoS for the d-wave Bogoliubov-de 
Gennes kicked rotator with decreasing disorder range as 
defined in Eqs. @, @ and @, with L x x L y = 256 x 256, 
A = 0.4, E F = 2tt 2 /5, K = 8., and N H = 1, 3, 5, 7, 17, and 
27 (from top to bottom). 

(black circles) which is in qualitative agreement with the 
prediction p(E) oc E~ l \ \iiE\- 3 of Ref. [13. 

Having established the validity of semiclassical predic- 
tions in the quantum chaotic regime, we next decrease 
the range of the disorder and enter the quantum disor- 
dered regime. We accomplish this via the inclusion of 
higher harmonics to the kicking potential, and replace 
Eq. © by 

KI N " 

^ K = ~\F2 — cos [Zx] cos [my] a z . (6) 

The typical impurity size decreases as C oc N^ 1 . Fig. 3 
shows the disappearance of the low-energy peak in the 
DoS as Njj increases. For the set of parameter con- 
sidered, once Nh — 27 is reached, the DoS vanishes at 
E = 0. Note that the resolution used in Fig. 3 does not 
allow to see the opening of the RMT gap below E* . A 
more precise look at the DoS for Njj — 27 is provided on 
the inset to Fig. 2 (empty diamonds). The data clearly 
indicate the expected RMT linear suppression of the DoS. 

Our results thus clarify the competition between RMT 
and semiclassics [ljj. The next step is to investigate 
the transport properties and to explore the parametric 
dependence of E* . Work along those lines is in progress 
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